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Abstract 

Using a large sample of photoproduced charm mesons from the FOCUS 

experiment at Fermilab (FNAL-E831), we observe the decay D° — ► K + ir~ 

with a signal yield of 149 ±31 events compared to a similarly cut sample 

consisting of 36 760 ±195 D° — > K~ir + events. We use the observed ratio of 

D°^K+it- to D°^K-tt+ (0.404 ±0.085 ±0.025)% to obtain a relationship 

between the D° mixing and doubly Cabibbo suppressed decay parameters. 
13.25.Ft, 12.15.Ff, 14.40.Lb 
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The decay D° — > K + 7r [[jj may occur either via a doubly Cabibbo suppressed (DCS) 



decay or through the mixing of the D° into D° followed by the Cabibbo favored (CF) decay 



D° — >K + 7i~ . The naive expectation for the ratio of DCS to CF branching fractions, -Rdcs ; 
is tan 4 #c" — 0.25%, which may be modified by final state interactions. Contributions from 
nonperturbative long range interactions make exact calculations of Standard Model charm 
mixing difficult, but the rate is expected to be less than 10~ 3 |2]J|. On the other hand, 
new physics may enhance mixing Q. Since Standard Model charm sector CP violation is 
expected to be small, and current searches report negative results f|, we ignore CP violation 
in this study. Conversely, recent charm sector mixing searches hint at an effect with a rate 
of order 10~ 4 |||7[]; hence possible mixing effects must be considered. 

Four groups |7 HT0|1 have studied the decay D° — > K + n~ and measured its branching ratio 
with respect to D° — > K~it + , but only the most recent result [|7| by CLEO II.V is statistically 
significant. In addition, there is a variation of a factor of five among the measurements with 
the most recent yielding the most events and the lowest value. We present a high statistics 
measurement of the branching ratio with more events than any previous experiment and 
different systematic uncertainties. 

The data were collected by the FOCUS Collaboration during the 1996-97 Fermilab fixed 
target run in the Wideband Photon beamline using an upgraded version of the E687 spec- 
trometer |TI|. Charm particles are produced in the interaction of high energy photons 
((E) ~ 180 GeV) with a segmented BeO target. In the target region, charged particles are 
tracked by 16 layers of silicon microstrip detectors which provide excellent vertex resolution. 
The momentum of the charged particles is determined by measuring their deflection in two 
oppositely polarized, large aperture dipole magnets with five stations of multiwire propor- 
tional chambers. Particle identification is determined by three multicell threshold Cerenkov 
detectors, electromagnetic calorimeters, and muon counters. 
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I. MEASUREMENT METHOD 



To minimize systematic effects we apply the same selection algorithm to both D° — > K + n~ 
and the normalizing mode, D° — > K~tt + . We use the sign of the soft tt (tt) from the decay 
D* + — > D°ff + to tag the production flavor of the D meson. For the normalizing mode the 
charge of the tt is the same as the charge of the D° daughter ir (right sign, RS), but is 
opposite for the DCS or mixing mode (wrong sign, WS). 

A D° candidate consists of a pair of oppositely charged tracks that form a decay vertex 
and have a Kir invariant mass between 1.7 and 2.1 GeV/c 2 . The D° candidate is used as a 
seed to locate a production vertex consisting of at least two charged tracks in addition to 
the D°. The production vertex is required to lie within la of the nearest target material and 
be separated from the decay vertex by at least five times the separation error {Ijai > 5). 
Each vertex must have a confidence level (CL) greater than 1% and the decay vertex tracks 
are required to be inconsistent with originating in the production vertex. 

Highly asymmetrical Kit pairs reconstructing with the D° mass are more likely to be 
background than signal. A cut is made on asymmetry, A = \pk~P-k I / \pk+Ptt I , to reject these 
candidates. The A cut point is lowered linearly as the D° momentum decreases to achieve 
the best rejection of background. 

For each charged track, the Cerenkov particle identification algorithm generates a set of 
X 2 -like variables, Wi = — 2 log(likelihood), where i ranges over the e, tt, K, and p hypotheses. 
We define W min as the smallest W of the four hypotheses and we require Wi — W min < 4, 
where % refers to either the K or 7r hypothesis. The D° daughters must also satisfy the 
slightly stronger Kn separation criteria of AWk = W^ — Wk > 0.5 for the K and AW^ = 
W K -W n > -2 for the tt. 

Events with the decay D° — > K~ir + where the K has been misidentified as a ir and the 
7r has been misidentified as a K, produce false WS candidates. These doubly misidentified 
events form a broad peak in the K + n~ mass distribution centered on the D° mass. When 
a real n tag is present, a peak indistinguishable from the real WS signal appears in the 
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D* — D mass difference. We treat this double misidentification background by imposing a 
hard Cerenkov cut on the sum AW 7 ^ + AW n > 8, when the invariant mass of the Ktt pair, 
with the K and tt particle hypotheses swapped, is within 4a of the D° mass. From Monte 
Carlo (MC) studies of this cut, we expect double misidentification feed-through in the WS 
signal to be consistent with zero, and we determine a 90% CL upper limit of 5% of the 
observed WS yield. 

All tracks assigned to the production vertex are considered as potential tt candidates. 
The tt candidate must satisfy W n —W m i n <4 and be inconsistent with being an electron from 
a 7 conversion where the 7 comes from a D°* decay. This is achieved using information from 
the Cerenkov system, electromagnetic calorimeters, and silicon microstrip detectors. 

Significant fake contributions to the WS yield arise from the decay D* + — >D°tt + in which 
the fr was correctly reconstructed but the D° was not. In the inset of Fig. [l|a we show the 
combined contributions from MC generated decays D° — * K + K~ and D° — > tt + tt~ , which 
are reconstructed as Ktt. Both the KK and tttt reflection peaks (below and above the D° 
mass respectively) have tails which extend into the D° mass region. In the inset of Fig. |T]b 
we plot the K + tt~ mass contributions from MC events of all known D° modes, except two- 
body final states. The flat background in the D° signal region is composed primarily of 
partially reconstructed and doubly misidentified D° decays to K~tt + tt° and K~t r v. The 
mass difference plots (Fig. |l|a and b) for reflected events in the D° signal region show peaked 
backgrounds in the D* signal region. Using tighter particle identification to eliminate these 
backgrounds rejects about one third of our real signal events. To avoid this we divide the 
RS and WS samples into 1 MeV/c 2 wide bins in mass difference from 139 to 179 MeV/c 2 
and plot the Ktt mass for each bin (a typical mass plot is shown in Fig. Q). The Kn mass 
distribution is then fit in each bin as follows: the structured reflections, KK and irn, are 
fit using line shapes obtained with the MC, the unstructured background is fit by a degree 
two polynomial, and the D° signal is fit with a Gaussian. By fitting in this way the real 
D° — > Kit decays are isolated from the correlated tt backgrounds. 

The D° yields from the 80 Ktt fits (40 each RS and WS) are plotted versus mass differ- 
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ence. The two composite mass difference distributions, shown in Fig. ^ are fit for the WS 
to RS branching ratio (-Rws)- The background is fit using the following function: 

/(Am) = a[(Am - m w ) 1/2 + /3(Am - m w ) 3/2 ], (1) 

where a and (3 are fit parameters and separate parameters (ckrs, Pus, a ws and /3ws ) are 
used to fit the RS and WS distributions. The shape of the high statistics RS signal is used 
to fit the WS signal. In the WS D* signal region the full fit function is the sum of the WS 
background parameterization and the scaled RS signal, where the signal scale factor is -Rws- 
Modeling the WS signal in this way avoids complications arising from parameterizing the 
non-Gaussian signal. We obtain .Rws = (0.404 ± 0.085)% and find 36 760±195 events above 
background in the RS signal region corresponding to a WS yield of 149 ±31 events. 

Several sources of possible systematic errors were investigated, taking care to avoid tests 
correlated to possible mixing effects. Since the WS and RS modes are kinematically identical, 
we expect the systematic effects due to spectrometer acceptance and analysis cuts to cancel. 

The effectiveness of the double misidentification cut was tested by watching the stability 
of i?ws as the Cerenkov cut is varied and the swapped-hypothesis mass window is widened. 
We found no evidence of a systematic error due to doubly misidentified decays. 

We investigated various methods of fitting the mass difference background such as us- 
ing a different parameterization, and constraining the RS and WS shapes to be the same 
(Pks — Pws)- We looked for sensitivity to the MC reflection shapes by shifting the reflection 
distributions by ±2 MeV/c 2 , and we searched for a systematic dependence on the mass 
difference binning by shifting the bin centers and changing the bin widths. All such tests 
returned values of -Rws consistent with the baseline measurement. 

To estimate the systematic error, measurements of -Rws were made with 136 different 
combinations of fit conditions and cut variations. Each measurement is assumed to be 
equally likely and we take the statistical variance of the measurements to be the systematic 
error on -Rws- We obtain a systematic error of 0.025%. 
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II. i?ocs IN THE PRESENCE OF MIXING 



The time dependent rate for WS decays relative to the CF branching fraction is 

R(t) = Rvcs + y't + {x ' 2+ / 2) t2 e ~*> ( 2 ) 

where t is in units of the D° lifetime. We use the strong phase (5) rotated convention of 
CLEO f7| , with y' = y cos 5 — x sin 5 and x' = x cos 5 + y sin 5 where x and y are the usual 
mixing parameters, x = Am/T and y = AT/2T. 

Using a large MC sample (10 times the data) of D° — > K~it + decays generated with a 
pure exponential lifetime of 413 fs |12 |, we can calculate the expected number of WS events 



by reweighting each accepted MC event with a weight given by 



Wi = ^ [ Jfccs + M^s y% + ^±A 2 1 , (3) 



N MC V — v 1 4 

where ti is the generated proper time for event i, and A^j a ta an d ^Vmc are the number of 
accepted RS events in the data and MC. Summing Eq. ^| over all accepted MC events and 
dividing by iV data we obtain 

Rws = Rdcs + V^r^ y'(t) + ^^(t 2 ). (4) 

The averages (t) and (t 2 ) are measured from the generated lifetime of the MC events accepted 
in the analysis. We find (t) = 1.578 ±0.008 and (t 2 ) = 3.61 ±0.03, where the errors are 
systematic, determined by comparing the reconstructed MC lifetime averages to the averages 
measured in data. Using our measured value of i?ws we obtain an expression for Rdcs as a 
function of the mixing parameters x' and y' . 

In Fig. [| we plot -Rdcs as a function of y' for two values of x' that cover the CLEO 
95% CL of \x'\ < 0.028. For comparison, the mixing measurements of CLEO and FOCUS 
are included in Fig. |j. The CLEO result comes from a direct measurement of -Rdcs? x ' an d 
y' . The FOCUS band represents a measurement of y from lifetime differences between CP 
even and CP mixed final states. The FOCUS y measurement can only be directly compared 
to the other results in the case of 5 = 0. For a model dependent comparison of the CLEO 



and FOCUS direct mixing results see |13 



III. CONCLUSION 



We observe a signal in the decay channel D° — > K + ir~ and measure its branching ratio 
relative to D° — > K~tt + to be (0.404 ±0.085 ±0.025)%. If charm sector mixing is significant, 
the doubly Cabibbo suppressed component of the branching ratio can be determined from 
the measured ratio by using the relation expressed in Eq. [| and plotted in Fig. [|. If charm 
mixing is sufficiently small, the doubly Cabibbo suppressed branching ratio is simply equal 
to the measured ratio. For comparison, Table [I] lists the existing measurements of this 
branching ratio, made with the assumption of no mixing. 
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FIGURES 

FIG. 1. Monte Carlo studies of background contamination from (a) D° — > K + K~ and 
—>7r + 7r , and (b) all known multibody D° decay modes. 

FIG. 2. A sample Kir mass fit. 

FIG. 3. (a) The RS mass difference distribution, with the inset showing a close-up of the RS 
background fit and signal region, (b) The WS mass difference distribution with the signal and 
background fit contributions shown. 

FIG. 4. Rbcs plotted as a function of y'. Contours are given for two values of x' covering the 
95% CL of the CLEO II. V result. 
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TABLES 

TABLE I. Measurements of i?DCS with the assumption of no charm mixing and no CP viola- 
tion. 



Experiment -Rdcs (%) n ° Mixing Events 
CLEO g 0.77 ±0.25 ±0.25 19.1 
E791 g 0-68lo;i ± 0.07 34 
Aleph || l-77±g:|g ± 0.31 21.3 

CLEO II.V 0.332±g;g|| ± 0.040 44.8 

This Study 0.404 ± 0.085 ± 0.025 149 
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